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1. Introduction 

The Feynman-Kac formula P provides a powerful connection between the imaginary 
time evolution of quantum systems in the continuum and probabilistic expectations over 
Wiener classical trajectories. Its extension to the case of multicomponent wavefunctions 
requires the introduction of a further expectation over stochastic Poisson processes 
[2]. The role of Poisson processes is central in the probabilistic representation of 
Berezin integrals over anticommuting variables and, in general, of the time evolution 
of discrete systems In the simplified formulation it is shown that the real 
time or the imaginary time dynamics of systems described by a finite Hamiltonian 
matrix, representing bosonic or fermionic degrees of freedom, is expressed in terms of 
the evolution of a proper collection of independent Poisson processes. For a lattice 
system the Poisson processes are associated to the links of the lattice and their jump 
rates can be arbitrary. In [Jj it is demonstrated that when the rates of the Poisson 
processes are chosen equal to the hopping coefficients of the system the probabilistic 
representation leads to an optimal algorithm which coincides with the Green Function 
Quantum Monte Carlo (QMC) method E| in the limit when the latter becomes 
exact jH]. 

In this paper we exploit the probabilistic representation [3] to derive analytical 
expressions for the matrix elements of the evolution operator in the long time limit. 
Our approach is based on a series expansion of the probabilistic expectation in terms of 
conditional expectations with a fixed number N of jumps of the Poisson processes. This 
resembles the expansion of the grand canonical partition function in statistical mechanics 
in terms of canonical averages with a fixed number of particles. By integrating out the 
N stochastic jump times, we show that the conditional expectations become averages 
of functions which depend only on the multiplicities Ny and of the values assumed 
by the potential and hopping energies, respectively, of the configurations visited by 
the system. According to a central limit theorem, at large values of N, i.e. long 
times, the rescaled multiplicities Ny/y/N and Na/VN become Gaussian distributed 
and the corresponding averages can be evaluated analytically. The parameters of the 
Gaussian probability density do not depend on the Hamiltonian parameters and are 
easily determined statistically. Finally, the series of the conditional expectations can be 
re-summed with a saddle-point method. As a remarkable result, we find that the ground 
state energy is semi-analytically determined by solving a simple scalar equation. Once 
this equation is solved the quantum expectation in the ground state for other operators 
can be determined analytically by using the Hellman-Feynman theorem. The result is 
valid for boson or fermion systems. As regards fermions the well known sign problem 
[HI EE El EEU U21 EE] is avoided by introducing an approximation related to the exact 
counting of the positive and negative contributions in the noninteracting case. 

The paper is organized as follows. In Section |21 we review the probabilistic 
representation of quantum dynamics in the case of a generalized Hubbard Hamiltonian. 
In Section El we decompose the expectation in canonical averages of weights which are 
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calculated analytically. The canonical averages are evaluated at long times via a central 
limit theorem in Section 0J The equation for the ground-state energy is discussed in 
Section 14.11 for hard-core bosons and in Section 14.21 for fermions. In Section 03 we show 
how to calculate ground-state correlation functions within our approach. Finally in 
Section IH1 we show some example cases and compare with the results of exact numerical 
calculations. General features of our approach are summarized and discussed in Section 

El 

2. Exact probabilistic representation of lattice dynamics 

We illustrate our approach in the case of imaginary time dynamics for a system of 
hard-core bosons or fermions described by a generalized Hubbard Hamiltonian 

H = - Yl viAvCj* + ^ c h c iAi c n + (i) 

i^jeAcr=TJ. i£A i£A a=U 

where A C Z d is a finite d-dimensional lattice with |A| sites and q ct the commuting 
or anticommuting destruction operators at site i and spin index o with the property 
cf CT = 0. We are interested in evaluating the matrix elements {n\e~ Ht \no) where 
n = (nif, nil, ■ ■ ■ > n AIT' n \Ml) are ^ e lattice occupation numbers taking the values 
or I. The total number of particles is N a = J2 ieA n ia ^ or a = tl- 111 the following we 
shall use the mod 2 addition n®n' = (n + n') mod 2. 

Let r be the set of system links, i.e. the unordered pairs with i, j e A such 
that rjij 0. For simplicity, we will start by assuming rjij = e if i and j are first 
neighbors and rjij = otherwise. We will also assume 7« = 7 and <5j = 0. We shall call 
such a model the first neighbor uniform (FNU) model. For a <i-dimensional lattice the 
number of links per spin component is |T| = d\A\. Let us introduce 

\ja{ri) = (n © l ia © \ ia \c\ a c ia + c) a c ia \n) (2) 
V(n) = (n\H\n), (3) 

where = (0, . . . , 0, l^, 0, . . . , 0), and let {N^ a }, G T, be a family of 2|T| 

independent Poisson processes with jump rate p. At each jump of the process N^ a , 
if Xij a 7^ a particle moves from site i to site j or vice versa, while the lattice 
configuration n remains unchanged if Ajj CT = 0. The total number of jumps at time t is 
j)ei>=u Njjer- By ordering the jumps according to the times s k , k — 1, . . . , N t , 
at which they take place in the interval [0, t), we define a trajectory as the Markov chain 
rii, n 2 , . . . , n Nt generated from the initial configuration n . Let us call Ai, A 2 , . . . , Ajv t 
and Vi, V2 . . . , Viv t the values of the matrix elements (J2J and (jHJ) occurring along the 
trajectory. As proved in jl], the following representation holds 

(n\e- Ht \n ) = E(5 n , nN M t ), (4) 

where the stochastic functional .M* is defined by 

M t = e 2\r\ P t j j-j £ Afce -v fc _ 1 (. fc - Sfc _ 1 ) j e -v Nt (t-s Nt ) (5) 
\fc=i p ) 
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if N t > and M l = e mpt e -v t if N t = 0. Here V = V{n ) and s = 0. 

Several quantities can be obtained from the matrix elements (j3J). The ground-state 
energy is given by 

Eo = Um -gf<^;>°> = , im =^i. (6) 

t^oo 22n\ n \ e ~ \ n o) E(M*) 
3. Canonical decomposition of expectation 

To evaluate © we decompose the expectation E(_M') as a series of conditional 
expectations with a fixed number of jumps (canonical averages) 



E(M t ) = J2 E (M*\N t = N) 



N=0 

OO 



= E E (7) 

N=oren N 

where — ^Jv( n o) is the set of trajectories with N jumps branching from the initial 
configuration no and 

cM _ \(r) \(r) x (r) /n\ 



W^ r) (t) = eM dsi / ds 2 ... dsNe -Vos 1 -vr(s 2 s l) -...-v^(t-s N )_ (Q) 

The weights (jHJ) are obtained on multiplying (jHJ by the infinitesimal probability 
e _2 ' r ' pt p 7V (isi(is2 ... (is at to have AT jumps and integrating over the jump times. 

A link ij with spin a is called active if \j a ^ 0. From (jHJ) it is clear that only 
trajectories formed by a sequence of active links contribute to ((2)). Hereafter we restrict 
VLn to be the set of these effective trajectories with A" jumps. The sum over the set Qn in 
(JJJ) can be rewritten as an average, (•), over the trajectories with A" jumps generated by 
extracting with uniform probability one of the active links available at the configurations 
no, rii, . . . , njv-i- If A-k = Y,n j)er a=u \^ijv( n k)\ is the number of active links in the 
configuration rik, the probability associated to the trajectory r is p^ = Ylk=o V^aT^ 
and we have 



N-l 

y: spwP( t ) = £ ppspwp® n 4 r) 

N-l 

(S N W N (t) 11 A k ) (10) 



fc=0 

rJV-l 



Note that Ylk=o = YI.A A Na depends only on the multiplicities Na of the values A 
assumed by the number of active links; these multiplicities are normalized to N, i.e. 
YIa Na — N . For the FNU model the possible values of A depend on the number of 
particles and we have the bound A < min(2<i(A r | + AjJ, 2|T|). 
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For a generic trajectory the weights (jHJ) satisfy the recursive differential equation 

= w N -i(t) - v N w N (t), (ii) 

where W_i(t) = 0. In terms of the Laplace transform Wn(z) = dte~ zt Wjsiit) we get 

^)^ N fl~^Uj^- (») 

k=0 V K ' 

In (112}) we see that the weights depend only on the multiplicities Ny of the values V 
assumed by the potential; these multiplicities are normalized to N + 1, i.e. ^y Ny = 
N + 1. For the model (JTJ) the possible values assumed by V are V = 0, 7, 27, . . . , N p j, 
where N p = min(iVj, iVj). 

For an assigned set of multiplicities of the potential, the antitransform of ()12j) can 
be evaluated with the residue method. In this way one obtains an exact recursive 
expression of WW(t) which, however, for large N must be evaluated numerically with 
multi-precision algebra. On the other hand, for N large a complex saddle-point method 
can be used which provides the following asymptotically exact explicit expression 

p x t-Z v N v log[(x +y)/e] 

W N (t) = ^ , (13) 

where x is the solution of the equation 

T_ J ^ = t. (14) 

Note that in the case 7 = 0, Eq. (fT3J) reduces to Stirling's approximation of the exact 
value Wjv(t) = e N t N /N\. 



4. Canonical averages via a central limit theorem 

As is evident from the explicit expression given in the case 7 = 0, the weights Wjv(£) 
have a maximum at some N which increases by increasing t. This remains true also 
in the general shown in the following. Therefore in the long time limit the 

most important contributions to the expansion (jJJ) of the expectation E(A^') come from 
larger and larger values of N. In this section we will evaluate the canonical averages (jl(Jj) 
analytically for iV large by using the asymptotic behavior of the stochastic variables Ny 
and iV^. In this limit we will not distinguish the different normalizations, N + 1 and 
N, of Ny and iV^, respectively. In the following, we will indicate with my and the 
number of different values assumed by the variables V and A, respectively. For clarity, 
we consider separately the hard-core boson and fermion cases. 



4-1. Hard-core bosons 

In this case we have Sn = 1 and the canonical averages ()1U|) are averages of a function 
which depends only on the multiplicities, Ny and Na (besides a parametric dependence 
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on time). In terms of the corresponding frequencies, vy = Ny/N and va — Na/N, which 
for N large become continuously distributed in the range [0, 1] with the constraints 

X> = 5>a = i, as) 

V A 

Eq. (|T0|) can be rewritten as 

N-l „ 

(W N (t)Y[A k )= dvV N (v)g N (t;v) (16) 

fc=0 ^ 

g N {t- v) = n = , (17) 



where v and u are vectors with m = my + components defined as v T = 
(. . . Vy v a . . .) and u T = (. . .— log[(x + V r )/e] log A . . .), respectively. Note 

that u depends on u through x = Xq{u). 

The probability density Pjv(k') is given by the fraction of trajectories branching 
from the initial configuration no and having after iV jumps multiplicities Ny = uyN 
and Na = vaN . For N large, it can be approximated in the following way. We rewrite 
the multiplicities as N v = ^f=o XvK^fc) an d n a = SfcLi XA{n k -i), where xv(n) = 1 
if V(n) = V and Xv( n ) — otherwise, and similarly for xa- Since the configurations 
rtj, form a Markov chain with finite state space, a central limit theorem applies to each 
rescaled sum Ny/y/N or Na/VN [E]. However, due to the constraints (|15j). the joint 
probability for these m rescaled sums is not Gaussian. Given an arbitrary set of my — 1 
VMike components and — 1 A-like components, the joint probability density is the 
product of a Gaussian density for this set of variables and two delta functions which 
take into account the constraints (|15|). For the frequencies u, therefore, we have 



V A 

where Tn^) is the normal density defined in terms of the vector v having the m — 2 
chosen components of v 



Here V and are the (m — 2)-component subvector and submatrix, respectively, of 

the mean value, 17, and the covariance matrix, SiV" 1 , of v. As discussed in Section |BJ 
the quantities V and S are easily measured by sampling over trajectories with a large 
number of jumps. 

By using ()18)) . the m-dimensional integral which appears in can be performed 
by the saddle-point method. Note that this integration method is asymptotically exact 
for iV large. Due to the constraints ([T5|). V satisfies the property J2v = Y1a^ a = 
while rows and columns of the VV, VA, AV, and AA blocks of S are normalized to 
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zero. By using these properties, in terms of v and S we get 



(W N (t)Y[A k ) 



k=0 



V (xo+V) 2 



where v sp is the saddle-point frequency defined by the equation 



S P — 77 



v + Su v 



(20) 



(21) 



In order to evaluate the expectation value E(.M*) we need to sum the series (JIJ). 
This can be done with a further saddle-point integration. According to Eq. (j2U|) the 
terms of this series are exponentially peaked at N = N sp , where N sp satisfies 



{t-N [(17, V ) + (En, v)]} j£ + (17, u) + i (Su, w 



= 0, (22) 
0...). We observe that, according to Eqs. (fl"4"|) 



with v T = (. . . (x + V")" 1 
and (|21|L the term {t — iV [(17, t>) + (Sw, t))]} vanishes for z/ = v sp so that the above 
condition reduces to 

, 1 



v, u) 



Sw,u 



= 0. (23) 

j/=j/sp j 7V=A rs P 

Equation (J23|) is a time independent equation which determines Xo\ v = v s P N=N s P as a 
function of 17 and X. According to Eq. (|14jl. this means that for = i^ sp , the quantity 
increases linearly with time so that Xq\ 1/=1/SP j^_^ sp becomes independent of time. 
In conclusion, a saddle-point integration with respect to iV of the series (J2J) provides 

X()t 

(24) 



E(M') 



By taking the time derivative of this expectation and using JBJ), we obtain that the 
ground-state energy of the hard-core boson system is 



E 



OB 



Equation (j22J) is, therefore, the equation for the ground-state energy. It defines E 0B in 
terms of V and S and explicitly reads 

-E 0B + V s 







-J2^vy> log f - 

V,V' ^ 

1/4 V 



^z7 A log(A) 



0-B 



V 



Eob + V 



log 



0-B 



V 



log (A) 



i^S^log(A) log (A' 



(26) 



A,A> 



In the case 7 = the ground-state energy E^ can be solved explicitly and one has 



Kb = -eexp 



v A log(A) + - ^AA> log(A) \og(A' 



A. A' 



(27) 



Ground state of many-body lattice systems: an analytical probabilistic approach 



8 



Note that Eq. (J27|) is a nontrivial formula for the ground state of a system of bosons 
interacting via a hard-core potential. With the above expression Eq. ()26|) can be written 
more compactly as 

io g (^) = x>io g (=^±r) 

_V Evy , log (=^±z) log (z^±z:) 



^v, A lo g ( - EoB + V )lo g (A). (21 

V,A V e / 



By using the bounds E^ < E b < 0, the scalar Eq. (J2%|) can be easily solved with the 
bisection method. 

The generalization of the above results to Hubbard Hamiltonians (0) with arbitrary 
parameters r], 7, d is straightforward. Equations (j20H25|) remain formally unchanged, 
however the vectors u, u, v and the covariance matrix S are modified in order to take 
into account all the possible values of the generalized potential V corresponding to the 
operator J2i & A li c l\ c ii c li c ii + J2ieA ^v=n ^^ c ^ c ia and all the possible values of the 
generalized kinetic quantities T = Arj/e, where now e is a unity of energy and rj/e is 
the dimensionless hopping value corresponding to the current jumping link. In fact the 
generalization of Eq. fTT)|) consists in replacing n^To 1 with Ylk=o(^kVk/ e)- Explicitly, 
now the vectors V and u are 

V T = (. . .V v V T ...) 

u T =(...- log[(x + V)/e] . .logT. . .) (29) 
and the generalized ground-state energy corresponding to 7 = 5 = reads 



e ob = -eexp 



^ v T log(T) + - ^2 ^T.r log(T) log(T') 

T T T l 



(30) 



Similar considerations hold in the case of Hamiltonians with arbitrary potential 
operators. 



4-2. Fermions 

In this case we have Sn = ±1 and the canonical averages (jlOj) are averages of a function 
which depends not only on the multiplicities, Ny and AOi but also on AL, the sign 
multiplicity related to by = (— l)^ - . The approach developed in the case of 
hard-core bosons can be extended also to fermions by including this further multiplicity 
iV_. We will report on this procedure elsewhere. Here, we introduce an approximation 
which allows to reduce the calculation of the fermion ground-state energy to that of an 
effectively modified hard-core boson system. This approximation is motivated by the 
observation that the correlations between N_ and Ny are smaller than those between 
N- and N A . 
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Let us consider again the FNU model. In order to evaluate (jlUj) for a fermion 
system we introduce the average weighted sign s^ after N jumps 



The quantity s^ is a function of the interaction strength. For 7 = it can be evaluated 
in the following way. Expanding X!n( n l e ~ m | n o) i n powers of t and comparing with 
the expansion ((7|), for 7 = and N large we get in the case of hard-core bosons and 
fermions, respectively 

N ~ 1 N 

(i[A k ) = c OB (n )(-E^/e) (32) 



N ^ N 



(S N H A k ) = c 0F (n ) (-£$/e) , (33) 

k=0 

where Cos (no) and cof(no) are coefficients related to the initial configuration no and 

and E^p are the 7 = grou 
jumps for 7 = is then given by 



and E^p are the 7 = ground-state energies. The average weighted sign after N 



8 „ = **pL e »H B ®''®). (34) 
coB(n ) 

In the case of fermions the noninteracting energy, E^p , is known exactly while for hard- 
core bosons can be computed with Monte Carlo simulations or analytically as shown 
above. In general E^p/E^ < 1 so that sn vanishes exponentially for N large. 

Approximating s^ with its value (|34j) at 7 = removes effectively the negative 
signs in the expectation E(Ai t ). Therefore this can be evaluated as in the case of hard- 
core bosons with the same Gaussian probability density. In particular the saddle-point 
condition for iV = N sp now becomes 

,(o)" 



1 E 
(F, «) + _(Su,u)+ log 



= 0, (35) 



which is a time independent equation determining xq\ v _ vSP n=n*v m the fermion case. 
Finally, as in the hard-core boson case one finds that 

Eof = ~ x o\„ =l ,s P jj =N sp , (36) 

and Eq. (|55j) becomes equal to Eq. (|28j) with E 0B and E^ substituted by E 0F and 
E^p, respectively. However, in this case we do not have the analogue of Eq. (|27|) and 
E^p must be provided by an independent calculation, i.e. by exact diagonalization of 
the separable many-particle Hilbert space. 



5. Ground-state correlation functions 



The equation obtained in the previous Section for the determination of the ground-state 
energy depends on the parameters of the Hamiltonian only explicitly through the values 
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of the generalized potentials V and of the kinetic quantities T. In fact, the statistical 
moments V and "S are determined by the structure of the Hamiltonian not by the values 
of the Hamiltonian parameters. Therefore we are able to evaluate the derivatives of 
the ground-state energy with respect to any parameter £ of the Hamiltonian H(£). 
This allows the determination of arbitrary ground-state correlation functions via the 
Hellman-Feynman theorem 

<s.K)i*«» — £ oK) ' ( ' 

where Eq(£) is the ground-state energy of H(£). Hereafter, we assume a normalized 
ground state, (£ (£)|£ (0} = 1- 

Suppose that we want to evaluate the quantum expectation of an operator O in the 
ground state of the Hamiltonian H. We have two possibilities. 

i) The operator O is a term of the Hamiltonian itself, e.g. O = c\ a Cj a , with i ^ j, 
O = 4t c it c lj, c ij,' OT O = c\ a c ia if H is the generalized Hamiltonian ((TJ). In this case, by 
using Eq. ()37|) for hard-core bosons we have 



(E OB (ri,-jf,d)\cl a c ja \EoB(il,l,S)) = ~ 9 mj E 0B (rt, 7, S) 

jT 1 



sp 



Ex 



i + J (3* 



'V -E 0B (v,f,S)+V 

(^0B(r7,7,<5)|cl T c iT 4^^1^05(77,7,(5)) = d^E^-q, 7, <5) 

sp 

V ^ <9 1/ 

_ -goa(T7,-7,g)+V ' ' 

17 -E 0B ( vn ,S)+V 

(E 0B (r7,7, (5)14^1^05(77,7,(5)) = d Sia E 0B (rj,j, S) 

V ^ & T/ 



(39) 



(40) 

V Yv_ 

l^V -E 0B ( v ,~y,5)+V 

The expressions in the second lines of (j38H40j) have been obtained by using the derivatives 
of Eq. ([23)1. which for a generic parameter £ read 

- £ »v d i log ( ~ EoB + V ) + £ 4^ log (T) = 0, (41) 
v \ e J T 

where v sp is given by (|21|) and is determined once S _b(?7, 7, 5) has been solved. Similar 
expressions hold in the case of fermions by using the derivatives of Eq. (|35|) 

- £ ^ log ( ~ jg0F + M + £ ^ CO + d, log f §|) = 0. (42) 

ii) If the operator O is not a term of the Hamiltonian H we consider a new 
Hamiltonian i?(£) = H+^O and calculate the corresponding ground-state energy Eq(£). 
Note that, since the used probabilistic representation holds for any system described by 
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a finite Hamiltonian matrix [3] , the nature of the operator O is arbitrary. As an example 
we study the spin-spin structure factor 

S(q x , q y ) = m E e^-^+^^HEolSiSjlEo), (43) 
' ' i,jeA 

where Si = cLc^ — c|, Qj, and Xj and are the coordinates of the i-th lattice point. The 
quantum expectation of the operators SiSj in the ground state of the hard-core boson 
FNU model, (E 0B (e, ^ISiSjlEos^e, 7)), can be obtained by considering the Hamiltonians 

H{£ ij ) = H + Z ij S i S j , (44) 

where H represents the FNU model. For these Hamiltonians, the possible values of the 
potential are V = 7717 + k£ij, with m — 0, 1, . . . , N p and k = — 1, 0, 1, and we have 

sp 

V ^ (9t V 



1^=0 



V 



,(45) 



(E OB (en)\SiSj\EoB(e,-f)) = <K _ /•-'()/;(' • " • G;) 

N 'V -i?os(e,7,fe)+^ 

where E s(e, 7, is the ground-state energy of the Hamiltonian (jjjj) and is calculated 
as explained in Section l4~Tl 



6. Numerical results 



Now we apply the approach developed in the previous Sections to some example cases. 
In particular, we compare the ground-state energy obtained by Eqs. (J28j) and (J35|) with 
that from exact numerical calculations. 

In our approach the starting point is the evaluation of the statistical moments V 
and Yi. These are obtained by generating trajectories in the lattice configuration space 
and counting the multiplicities Ny and Na- The length of the trajectories is chosen to be 
sufficiently large for the asymptotic behavior to be established. The determination of V 
and I] with good statistical precision requires a number of trajectories which increases 
no more than linearly with the number of lattice sites |A|. Therefore the evaluation 
of these moments is feasible even for large systems. In the following applications the 
statistical errors associated to the measurement of V and S are negligible on the scales 
considered. 

In Fig. we show the behavior of Eqb for the hard-core boson FNU model as a 
function of the interaction strength 7 for several lattice systems. The solution of Eq. 

compares rather well with the results of exact diagonalizations or QMC simulations. 
There is a small systematic error which grows with increasing 7 and/or the system size 
and which is maximum at half density. This error is related to the Gaussian shape 
fll9jl assumed for the asymptotic probability density. In fact, the mentioned central 
limit theorem for Markov chains applies to the variables vy\f~N and va^/N while the 
function given by Eq. (JT7j) depends on vyN and vaN . This implies that the tails of 
the probability density give a finite contribution to the integral (|2L)j) . From the structure 
of gjsr, furthermore, it is evident that this error becomes large when the components of 
the vector u assume large values, i.e. for large values of 7 and/or large system sizes. 



Ground state of many-body lattice systems: an analytical probabilistic approach 12 



0| , . . , . , . , , . r 




4 8 12 

7 A 

Figure 1. Ground-state energy per particle for the hard-core boson FNU model vs 
the interaction strength 7. The results from Eq. I|28[l (solid lines) are compared 
with those from exact diagonalizations (+) and QMC simulations (x) for two- 
dimensional systems with periodic boundary conditions and different L x x L y sites 
and iVf = N± = N p particles. The statistical errors in the QMC data are negligible in 
this scale. 




4 8 12 

7 A 

Figure 2. As in Fig. Q]f or the fcrmion FNU model. In this case no QMC simulations 
are available. Exact diagonalization data for the 4x4 systems are taken from ^] 
(N p = 8) and PSE] (N p = 5). 

In Fig. |21 we show the behavior of E p evaluated according to Eq. (|35|) as a function 
of the interaction strength 7 for the same cases considered in Fig. Compared with 
the hard-core boson case we observe a further systematic error due to the approximation 
sjv(7) — sjv(O). Depending on the particle density of the system, this error adds to or 
subtracts from the systematic error due to the Gaussian tails of the probability density 
discussed in the case of hard-core bosons. 

We stress that, once V and S are known, our approach provides any other ground- 
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Figure 3. Spin-spin structure factor S(q x ,q y ) for the hard-core boson FNU model 
at 7 = in a lattice with 4x4 sites, iVf = N± = 8 particles, and periodic boundary 
conditions. For the same system, in the right plot we show the values of the ground- 
state expectation of the operators SiSj for dij = 0, 1, 2 as a function of the interaction 
strength 7 (solid lines) compared with Monte Carlo results (dots with error bars). 

state quantity analytically as a function of the Hamiltonian parameters. On the other 
hand, QMC methods require, due to the unavoidable branching or reconfiguration 
techniques jZj, different simulations for different values of the parameters. 

As an example of calculation of correlation functions in the ground state, in Fig. 
El we report the spin-spin structure factor S(q x , q y ) evaluated by using Eqs. ()43II45|) 
for the hard-core boson FNU model. The value S(tt, it) represents a maximum for 
S(q x , q y ) and for large values of |A| is related to the staggered magnetization m s through 
m s = a/ £(71-, 7r) / 1 A I . Note that in the 4x4 system with periodic boundary conditions 
considered in Fig. El we have only five different values for (Eo\SiSj \E ) corresponding to 
the five possible distances = \xi — Xj \ + |y$ —yj \ = 0, 1, 2, 3, 4. In Fig. El we also show 
the behavior of (Eo\SiSj\Eo) as a function of 7 for = 0, 1,2. These terms provide 
the most important contributions to S(ir, 7r). For small values of 7, the results compare 
rather well with Monte Carlo data. 

7. Conclusions 

By using saddle-point techniques and a central limit theorem, we have exploited an exact 
probabilistic representation of the quantum dynamics in a lattice to derive analytical 
approximations for the matrix elements of the evolution operator in the limit of long 
times. For both hard-core boson and fermion systems, this development yields to a 
simple scalar equation for the ground-state energy. This equation depends on the values 
of the generalized potentials V and of the kinetic quantities T, and on the statistical 
moments u and XI of their asymptotic multiplicities Ny and Nt- In turn, these moments 
depend only on the structure of the system Hamiltonian, not on the values of the 
Hamiltonian parameters. This implies that the statistical moments must be measured 
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una tantum for a given Hamiltonian structure and, once V and S are known, our 
approach provides the ground-state energy analytically as a function of the Hamiltonian 
parameters. 

In the long time limit, the saddle-point integrations used in our approach are 
asymptotically exact and the central limit theorem evoked applies rigorously to the 
rescaled multiplicities Ny / y/N and Nt/VN. However, since functions depending on Ny 
and Nt are involved, we have a small systematic error related to the finite contributions 
from the tails of the probability density. This systematic error could be reduced by 
a large deviation analysis. In fact, equations (|27)l and ([2~H]) suggest that the Gaussian 
approximation for the probability density corresponds to a second order truncation 
of a cumulant expansion. Anyway, the present Gaussian approach has the following 
relevant features: i) the equations derived for the ground-state energy and the ground- 
state correlation functions are particularly simple; ii) the corresponding results compare 
rather well to the exact ones in regions of physical interest. 

In this paper we have considered two-dimensional lattice models at imaginary times. 
Our approach, however, is valid in any dimension. Similar analytical expressions can be 
obtained also for the real time evolution. 
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